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Abstract 



Simple, self-similar, analytic solutions of 1 + 1 dimensional relativistic hydrody- 
namics are presented, generalizing Bjorken's solution to inhomogeneous rapidity 
q ! distribution. 
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Relativistic hydrodynamics has various applications, including the calculations 
of single-particle spectra and two-particle correlations in high energy heavy ion 
collisions. The applications of relativistic hydrodynamics in the field of heavy 
ion physics have been reviewed in [1] . More recently, there has been an increas- 
ing interest in studying RHIC and coming LHC experiments in the framework 
of relativistic hydrodynamics [2-4]. The hydrodynamical analysis can also be 
extended to the study of these processes on event-by-event basis [5,6]. How- 
ever, most works in hydrodynamics are numerical so not always transparent. 
In this sense, exact solutions would be useful, but are rarely found due to 
the highly non-linear nature of relativistic hydrodynamics. Actually, Landau's 
one-dimensional analytical solution of relativistic hydrodynamics [7] gave rise 
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to a new approach in high energy physics. The boost-invariant Bjorken so- 
lution [8], found more than 20 years later, is frequently utilized as the basis 
for estimations of initial energy densities in ultra-relativistic nucleus-nucleus 
collisions. Recently, Biro has found self-similar exact solutions of relativistic 
hydrodynamics for cylindrically expanding systems [9,10], however, his solu- 
tions are valid only when the pressure is independent of space and time, as e.g. 
in the case of a rehadronization phase transition in the middle of a relativistic 
heavy ion collision. 

Here we present an analytic approach, which goes back to the data-motivated 
exact analytic solution of non-relativistic hydrodynamics found by Zimanyi, 
Bondorf and Garpman (ZBG) in 1978 for low energy heavy ion collisions with 
spherical symmetry [11]. This solution has been extended to the case of el- 
liptic symmetry by Zimanyi and collaborators in [12]. In [13,14] a Gaussian 
parameterization has been introduced to describe the mass dependence of 
the effective temperature and the radius parameters of the two-particle Bose- 
Einstein correlation functions in high energy heavy ion collisions. Later it 
has been realized that this phenomenological parameterization of data corre- 
sponds to an exact, Gaussian solution of non-relativistic hydrodynamics with 
spherical symmetry [15]. The spherically symmetric self-similar solutions of 
non-relativistic hydrodynamics were obtained in a general manner in [16], 
that included an arbitrary scaling function for the temperature profile, and 
expressed the density distribution in terms of the temperature profile func- 
tion. The ZBG solution and the Gaussian solution of [15] are recovered from 
the general solution of [16] as special cases, corresponding to different scaling 
functions of the temperature profile. The Gaussian solution has been general- 
ized to ellipsoidal expansions in [17], that provides analytic insight into the 
physics of non-central heavy ion collisions [18]. 

Our approach corresponds to a relativistic generalization of these recently ob- 
tained analytic solutions [15,16,18] of non-relativistic fireball hydrodynamics 
to the case of relativistic flows, based on the success of the analytic approach 
to parameterize the single particle spectra and the two-particle Bose-Einstein 
correlations in high energy heavy ion physics in terms of hydrodynamically 
expanding sources [19]. 

In particular, we attempt here to solve the 1+1 dimensional relativistic hy- 
drodynamical problem, in trying to overcome two shortcomings of Bjorken's 
well-known solution. These two shortcomings of Bjorken's solutions are that 
i) it only describes ultra-relativistic limit so that the rapidity distribution is 
flat; ii) it contains no transverse flow. 

Here we present a new family of exact analytic solutions of relativistic hydro- 
dynamics in 1 + 1 dimension (time + longitudinal coordinate), that is able to 
describe arbitrary inhomogeneous rapidity distributions. The inclusion of the 
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transverse flow will be treated in a subsequent paper. 

An interesting aspect of this solution is that the shapes of the rapidity distri- 
bution dN/dy and temperature distribution are coupled in the way that the 
larger the rapidity density, the smaller the effective temperature. Choosing 
the effective temperature distribution T e ^{y) to be flat, we recover Bjorken's 
1+1 dimensional solution. 



2 The equations of relativistic hydrodynamics 



In order to get solutions that have simple non-relativistic limiting behavior, 
we consider a gas of one type of particles whose total number is conserved. 
We then solve the relativistic continuity and energy-momentum conservation 
equation: 



c^(nO = 0, (1) 
d^ v = 0. (2) 

Here n = n(x) is the number density, the four- velocity is denoted by = 
(x) = 7 (1, v), normalized to u^u^ = 1, and the energy-momentum tensor 
is denoted by T^ u . We assume perfect fluid, 

T» v = (e + p) u»u v - pg^, (3) 



where e stands for the relativistic energy density and p denotes the pressure. 

We close this set of relativistic hydrodynamical equations with the use of an 
equation of state. We assume an ideal gas that contains massive conserved 
quanta, 



e = mn + up, (4) 
p = nT. (5) 

This equation of state has two free parameters, m and k. The case of ultrarel- 
ativistic gas is recovered by choosing m = and k = 3. If the thermal motion 
of the gas is non-relativistic, one has m ^> T. Using this condition, the equa- 
tions of non-relativistic hydrodynamics are re-obtained from the above set of 
equations in the v 2 < 1 limiting case. The non-relativistic ideal gas equation 
of state corresponds to the case of k — 3/2. 

The energy-momentum conservation equation contains 4 independent equa- 
tions. These can be projected into a component parallel to w M and a component 
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orthogonal to w M , that yield the relativistic energy and Euler equations: 



u"V+(£+ P )a/=o, (6) 

u^d^p + (e + pj/^w, - d uP = 0. (7) 

The relativistic Euler equation contains only three independent components, 
by construction. 

Using the equation of state and the continuity equation, the energy equation 
can be rewritten as an equation for the temperature, 

+ -TduU» = 0. (8) 

K 

From now on we assume that the flow is homogeneous in the transverse di- 
rection corresponding to infinitely broad target and projectile in high energy 
heavy ion collisions. We reduce our equations to the 1 + 1 dimensional case, 
so the coordinates are x^ = (t,r z ), x^ = (t,—r z ) and the metric tensor is 
<T = <7^ = diag(l,-l). 

We solve 3 independent equations, the continuity, the temperature equation 
and the z component of the Euler equation (1,7,8). The equations, (4,5) close 
this system of equations in terms of 3 variables, n, T and v z . 



3 Self-similar solution 



We search for solutions that scale in the z direction. We introduce the time 
dependent scaling variable 

x = — (9) 



and assume that the longitudinal motion corresponds to a Hubble type of 
self-similar longitudinal expansion, 

v z (t,r z ) = ^r z , (10) 



where Z = dZ(t)/dt stands for the time derivative of the scale parameter Z. 
In a relativistic notation, this form is equivalent to 
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«" = (coshC,sinhC), (11) 
Z(t) 

tanhC = 777T r ^ (12) 
s Z(t) 

cosh C = —j= 1 , =7- (13) 
VI - Z 2 a; 

Using this ansatz, we find that the continuity equation is solved by the form 

"('^""f (14) 



where Q (x) is an arbitrary non-negative function of the scaling variable x and 
n and Z are normalization constants. We use the convention Z = Z(t ) and 
n = n(to, 0) which implies that Q{x = 0) = 1. The temperature equation, (8) 
is solved by the following form: 




T(t,r M )=T ^— — F(x). (15) 



The constants of normalization are chosen such that T = T(t , 0) and .F(O) = 
1. Here again, we find that the solution is independent of the form of the 
function J-{x), provided that J-{x) > 0. 

Using the ansatz for the flow profile and the solution for the density and 
the temperature, the relativistic Euler equation reduces to a complicated non- 
linear equation that contains Z, Z and Z and x. Taking this equation at x = 
we express Z as a function of Z and Z. Substituting this back to the Euler 
equation we obtain an equation for Z, Z and x. In particular, for the m = 
case, Z cancels out and this reduces to a second order polynomial equation 
for Z 2 , which has only one positive root. The form of the solution in this case 
(m = 0) is 

Z\t) = F(x). (16) 



Observing that the function F depends only on the scaling variable x, while 
Z depends only on the time variable t, we conclude that the only solution of 
this equation should be a constant 

Z = Z . (17) 



Now we choose the origin of the time axis such that Z(t = 0) = without loss 
of generality. The solutions can be cast in a relatively simple form by intro- 
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ducing the longitudinal proper time r and the space-time rapidity r] defined 

as 



zi 

1, ft + r z 

^2 log U^ 



This implies that 



(18) 
(19) 



Z(t) = Z t, 

v z = -j- = tanh?7, 

C = v- 



(20) 
(21) 
(22) 



Thus the solution for the flow velocity field corresponds to the flow field of 
the Bjorken solution. However, in the Bjorken solution the temperature dis- 
tribution was independent of the r\ variable, while in our case the density and 
the temperature distributions can be both t] dependent, or in other words, our 
solutions are scale dependent. The scale is defined by the parameter Z , in the 
longitudinal direction. 

This special form of the solution for the flow velocity field implies that Z = 
0. This equation implies that there is no pressure gradient and there is no 
acceleration in this class of self-similar solutions. The Euler equation is reduced 
to the following requirement: 



(d, + -fd t ) 



1 - Zlxf +1 l^g(x)F(x) 



= 



(23) 



This equation is solved by the trivial Q{x)J r (x) = as well as by the non-trivial 
solution of 

Q{x)T{x) = (1 - Z 2 x)-^ K \ (24) 



which is indeed only a function of x as Zq is a constant of time. With this 
form, the Euler equation is satisfied. This solution implies that the scaling 
profile functions for the temperature and the density distribution are not in- 
dependent. As the constraint is given only for their product, one of them can 
be still chosen in an arbitrary manner. 

At this point it is worthwhile to introduce new forms of the scaling functions. 
Let us define 
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T(x)=F{x){l- Zlx) l/K , (25) 

V{x)=Q{x){l- Zlx). (26) 

Then the constraint Eq.(24) can be cast to the simplest form of 

V{x)T{x) = 1. (27) 



Let us summarize our new family of solutions of the 1+1 dimensional rela- 
tivistic hydrodynamics by substituting the results in the density, temperature 
and pressure profiles. 



We obtain 



T 

v z = j = tanh??, (28) 
n(t,r z )=n -V(x), (30) 

T 

,f \ 1+1/k 

P(t,r z )=p m , (31) 



T 



v(x)' 



where po = n T . 



This implies that we have generated a new family of exact solutions of rel- 
ativistic hydrodynamics: a new hydrodynamical solution is assigned to each 
non- negative function T(x). It can be checked that the above solutions are 
valid also for massive particles, the form of the solution is independent of the 
value of the mass m. The form of solutions depends parametrically on k, that 
characterizes the equation of state. 



4 Analysis of the solutions 



We have obtained new solutions of the (1+1) dimensional relativistic hydrody- 
namical equations which describe a self-similar, streaming flow. The pressure 
and the flow profiles are the same as in the 1+1 dimensional Bjorken solution. 
In the case of V(x) = 1, we recover Bjorken's solution. In this limiting case, 
the pressure, the density and the temperature profiles depend only on the 
longitudinal proper time r. 
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In the general case, our solution contains a characteristic scale defining pa- 
rameter in the longitudinal direction, Z , and an arbitrary scaling function 
V(x). Thus we have an infinitely rich new family of solutions. Let us try to 
determine the physical meaning of the scaling function V(x). 

In order to do this we evaluate the single particle spectra corresponding to 
the new solutions. Here we neglect any possible dynamics in the transverse 
directions, as usual in case of applications of Bjorken's solution. The four- 
velocity field of our solutions thus becomes = (cosh 77, 0, 0, sinh??). The four- 
momentum of the observed particles with mass m is denoted by k^ 1 = (m t cosh 
y, k x , k y , m t sinh y). Let us assume that particles freeze out at a constant lon- 
gitudinal proper-time 77, for the sake of simplicity. This implies freeze-out 
at a constant pressure, but at a space-time rapidity dependent temperature 
and density, and makes it possible to continue the calculation analytically. 
The source function of locally thermalized relativistically flowing particles in 
a Boltzmann approximation can be written as 

S(x, k) = C(rj) m t cosh(r] — y) n(x) exp (—k^u^/T) 5{r — r/), (33) 

where C(i]) is an rj dependent normalization factor, given by the condition 
that / dk/E S(x, k) = n(x)5{r — tj), which implies that 

C(rj) = {47rm 2 T(r / ,r ? )K 2 [m/T(r / ,7 ? )]} _1 , (34) 

where K v {z) = / °° dzexp(— z cosht) cosh(z/£) is the modified Bessel function 
of the second kind. 

The single particle spectrum can be calculated from the emission function as 
E— = / rdrdr]S(x,k). (35) 

Substituting our family of new solutions, and using T(x) = 1/V(x), we obtain 
the following form 



S(x, k) = C(rj) m t cosh(^ — y) n(x)fB{x, k) 



f B (x,k) =exp 



m t cosh(?7 — y) 



1 / K tanh 2 f], 
^0 



(36) 

5(r-r f ). (37) 



We are interested in the coupling between the measurable rapidity distribution 
and the rapidity dependence of the effective temperature in the transverse 
directions as obtained from our new family of solutions. We assume that V(x) 
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is a slowly varying function, i.e. dlog V(x)/dx <C 1 in the region of interest. 
This assumption implies that the point of maximal emissivity is located at rj = 
y with correction terms of 0{d\ogV{x)/dx) The measurable single-particle 
spectra can be written as 



E ^ = 2C (v) n ot V (^^j Kdmt/Taiy)], (38) 
dN „AanhV\ 

-dy- =notoV {-^r)- m 

where 



Tes(y) 



V 



( tanh 2 



Hi 



1/k 



(40) 



Note that the V function is a free fit function that describes the measurable 
rapidity distribution, including characteristic scales of the size of Z . 

We see that the slope parameter for transverse mass distribution T e g is related 
to the rapidity distribution as 

t ( ^ t dN/dy(y = 0) 

T c M = T [y f j dN/dy . (41) 



Figures 1 and 2 illustrate the calculated behavior of the effective temperature 
distribution as a function of rapidity for a single Gaussian-like and a dou- 
ble Gaussian-like ansatz for the measurable rapidity distribution. In case of 
a homogeneous rapidity distribution, dN/dy = C we recover Bjorken's result 
that the effective temperature distribution is rapidity independent. This be- 
havior is expected to appear in high energy heavy ion collisions in the infinite 
bombarding energy limit. 



5 Summary 



We have found a new family of solutions of (1+1) dimensional relativistic 
hydrodynamics. This family solves the continuity equation and the conserva- 
tion of the energy - momentum tensor of a perfect fluid, assuming an ideal 
gas equation of state. The flow field coincides with that of Bjorken's solution. 
However, the shape of the measurable rapidity distribution, dN/ dy plays the 
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role of an arbitrary scaling function in our solution, and we obtain that the 
effective temperature of the transverse momentum distribution becomes ra- 
pidity dependent. Assuming that dN/dy is a slowly varying function of the 
rapidity y, we find that the effective temperature is proportional to the inverse 
of the rapidity distribution, T e s(y) oc (dN / dy)^ 1 . 

As compared to the well-known Bjorken's solution, we have solved one more 
equation, the continuity equation. We have considered equations of state that 
have two free parameters, the mass m and k = de/dp. Interestingly, these 
generalizations resulted in additional freedom in the solution. Our solution, 
similarly to Bjorken's case, describes scaling longitudinal flow and a pressure 
distribution that depends only on the longitudinal proper time. However, in 
our case, the pressure is a product of the local number density and the local 
temperature, hence one of these can be chosen in an arbitrary manner. In 
principle, we obtained that the measured single particle rapidity distribution 
can be arbitrary, and this measurable distribution plays the role of a scaling 
function of both the density and the inverse temperature in this family of 
new solutions of relativistic hydrodynamics. A generalization of the present 
one dimensional solution to 3 dimensional case can also be done in a similar 
manner to the nonrelativistic case [17]. 
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Fig. 1. Rapidity distribution dN/dy and effective temperature distribution T e ^{y) 
as a function of rapidity y, as obtained from a new family of solutions of 
(1+1) dimensional relativistic hydrodynamics. Here we use the scaling function 
V(x) = (1 — xp l ^\ using a scale parameter Z$ = tanh(4), no^o = 900 and 
7b(W r /) 1/K = 200 MeV > corresponding to a single maximum in the rapidity distri- 
bution dN/dy. The analytic expressions are given by eqs. (30, 32,39,40). 




Fig. 2. Same as Fig. 1 but utilizing a different form of the scaling function, 
V(x) = y/l + 1.6x 4 — 2.6x 8 , using a scale parameter Zq = 1, ra io = 800 and 
Tb(^o/ r /) 1 ^ K = 200 MeV, corresponding to a two-peaked rapidity distribution. 
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